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We investigate the perturbation theory of a fixed-norm, timehke Lorentz-violating 
vector field. After consistently quantizing the vector field to put constraints on its 
parameters, we compute the primordial spectra of perturbations generated by infla- 
tion in the presence of this vector field. We find that its perturbations are sourced 
by the perturbations of the inflaton; without the inflaton perturbation the vector 
field perturbations decay away leaving no primordial spectra of perturbations. Since 
the inflaton perturbation does not have a spin-1 component, the vector field generi- 
cally does not generate any spin-1 "vector-type" perturbations. Nevertheless, it will 
modify the amplitude of both the spin-0 "scalar-type" and spin-2 "tensor-type" per- 
turbation spectra, leading to violations of the inflationary consistency relationship. 



I. INTRODUCTION 

Lorentz invariance is a cornerstone of both the standard model of particle physics and 
general relativity, two very successful descriptions of the observable universe. The former 
describes at the quantum level basic processes involving particles and forces, while the latter 
describes at the classical level long-range gravitational forces. It is believed that both theories 
will merge to form a complete quantum theory of gravity at some very high cut-off energy 
scale like the Planck-scale 
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This cut-off scale marks the point where our old description of nature breaks down, and 
it is not inconceivable that one of the victims of this breakdown is Lorentz invariance. It 
is thus interestin g to test the robustness of this symmetry at the highest energy scales 

The early universe provides such a laboratory to us. In this paper, we study the effects 
of such a fixed-norm, constant timelike Lorentz- Violating (LV) vector field (introduced in 
12 . 13 . 1^ ) during inflation and its imprints on the Cosmic Microwave Background (CMB). 



In order to do this, we extend the background theory discussed in a companion paper |15| 
into the regime of linear perturbations, investigating both the consistency of its quantum 
field theory and its evolution in an inflationary universe. 

Furthermore, in the theory of inflation, very short distance quantum fluctuations are 
"blown-up" by the accelerated expansion, generating a spectrum of primordial super-Hubble 
size fluctuations. Thus it is reasonable to ask whether the quantum perturbations of this 
vector field may undergo a similar procedure to leave us with a classical spectrum of perturba- 
tions. Combined with the fact that vector perturbations can source both B-type and E-type 
polarization modes [3] , such a primordial spectrum of vector modes could conceivably leave 
strong signatures on the CMB. 

The question now becomes whether perturbations of this vector field will be frozen-out 
during the inflationary epoch, akin to that of the inflaton and gravitational wave perturba- 
tions. As we will see in Section flVI|) . the answer to this question is no. 

Nevertheless, while we are thwarted in our attempt to directly see this vector field, it does 
leave other more subtle effects on the CMB. We investigate these effects in Section ()VII|) . 
where we show that the vector field modifies both the spin-0 "scalar" and spin-2 "tensor" 
perturbation spectrum in possibly detectable ways. It follows that the so-called inflationary 



17| is also violated. 



consistency relation 

This paper is organized as follows. In Section (jTlJ, we state our theory and derive the 
basic equations of motion. Since the unperturbed field is a constant, we consider its classical 
perturbations in Section pilj) . We then quantize this theory in Section (jIVj) to show that the 
quantum field theory has positive definite Hamiltonian and its excitations are not tachyonic 
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for some suitable choices of the theory parameters. In Section (j^, we put the field into 
a Friedman-Robertson- Walker universe, and derive the basic linear order equations for the 
vector field. We put these equations to work in Section ()VI|) . solving the linear order equa- 
tions in a de-Sitter universe. In Section ()VII|) . we add an inflaton field and then compute 
the generated primordial perturbation spectra. Finally we conclude in Section (jVIIIIl . 



II. THE EQUATIONS OF MOTION 

We begin by deriving the basic equations of motion for the vector field theory. We will 
follow the notation of |l5|, where a more careful derivation is presented. In this paper we 
set ^ = c = 1; we call the gravitational constant that appear in the action G^, since as we 
have shown in it differs from the observed Newton's constant Gat. 

The action for the most general theory of a single fixed-length vector field whose deriva- 
tives are at most second order are given by Jacobson and Mattingly as 

where 

Cu = -PiW^u'W^u, - (32iV^,u^f - /SsV^u'W.u^ + X[u^u^ + m\ (2) 

Here the /3j's are dimensionless parameters of our theory, and A is a lagrange multiplier 
field; the vector field itself has a dimension of mass. We have dropped all terms higher than 
quadratic order in u. 

Although this theory is manifestly Lorentz- invariant, its always non- vanishing vacuum 
violates Lorentz symmetry by picking out a dynamical frame at every point in spacetime. 
One can say that it has spontaneously broken Lorentz symmetry since is non- 

vanishing, there is always a dynamically preferred frame picked out by the vector field at all 
points in spacetime. 

To see Lorentz- violating matter dynamics, in principle we need to couple this vector field 
to the matter action (see for example In this paper we are interested in its gravitational 

effects as a matter field. 
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The equation of motion for A enforces the fixed time-hke norm constraint 

g^uU^u" = -m^ (3) 

m being the squared norm of the vector field, while the equation of motion for the vector 
field itself is 

V^J^^ = Xu". (4) 

where the "current" tensor 

= -A V^M, - f32VXSii - /^aV.n^ (5) 
By multiplying Eqn. (@)) with on both sides, we obtain A = —m~'^Uu'V^J^'^. Using this 
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we find that the stress-energy tensor for this theory is 
Ta(3 = 2/3i(V«m^V/3Mm- V^n„V^M/3) 

-2m~'^u^Vi,J'"'uaUi3 + gapCu- (6) 

III. CLASSICAL PERTURBATIONS IN FLAT SPACE 

Before we study the behaviour of the vector field in a cosmological setting, we first consider 
the behaviour of this vector field in Minkowski space without gravity^. Since the unperturbed 
vector field is a constant, we investigate the evolution of its perturbations. In particular we 
want to put constraints on the parameters jSi, by appealing to the physical properties of its 
perturbations both classically and quantum mechanically. 

Our goal in this section is to find the classical solutions to the equations of motions of the 
vector perturbations. Along the way, we derive the Lagrangian for the vector perturbations, 
which we will use in Section pV|) when we consider the quantum version of this theory. 

In Minkowski space (with metric signature +2) without gravity, a solution to the equations 
of motion (jH)) and (jlj) is simply 

u^ = (m, 0,0,0) (7) 



^ This limit corresponds to the setting IGTrG^/Jim^ ^ 1, see Section (jVI D|l below. 
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and 

A = 0. (8) 
To investigate its dynamics we proceed to linear order by perturbing u and A 

— >u^' + v\ A — >X + SX. (9) 

We insert these expansions back into the Lagrangian Eqn. (^J, and expand to quadratic 
order 

Cu^Ca + 5iC + S2C (10) 

where 6iC and 62C contain terms of hnear and quadratic order in v respectively. Setting 
6iC = will recover the equations of motion for the background variables u, while variation 
of the quadratic order Lagrangian with respect to the perturbed variables will give us the 
first order perturbation equations of motion [2^. 

Since we are interested in the dynamics of the perturbations, we now focus on the 
quadratic order Lagrangian 

62C = 2 [-(3,dXd''vu - f32{dX? - Psd^v'd.v^ + 6X{u^v^)] . (11) 
where we have used the background result A = 0. The equation of motion for 6X is 

n% = 0' (12) 
and by recalling -u^ = (m, 0, 0, 0), we obtain 

= 0. (13) 

The vector field has no timelike perturbation; this follows from the fact that it is already 
constant in that direction by the fixed norm constraint and our choice of ansatz for its 
solution. 

Since SX is a lagrange multiplier, we can insert the result Eqn. ()13|1 back into the La- 
grangian Eqn. (|TT|l to obtain 



S2C = 2 



Piv' - (3,dyd'v, - I32{d,vy - mv'd.v' , (14) 
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where the Roman indices run from 1 to 3 and overdot denotes derivative with respect to time 
t. It is now perfectly reasonable to vary this action with respect to f * to obtain its equation 
of motion, since it is already in a respectable representation of the SO (3) group. However, 
this is a reducible representation, so under spatial rotations the components of the vector 
field will transform into each other. Thus it is more illuminating if we first decompose 
into its irreducible spin-0 and spin-1 components by 

yi = Qi^ + Ni (15) 

where V is a "scalar" potential and N'^ is the component transverse to d^V (i.e. diN"^ = 0). 

There are two reasons we do this. Firstly it extracts out the real dynamical degrees of 
freedom for the vector field, allowing us to quantize it separately. Secondly, we would like 
to consider its evolution in a cosmological setting, and decomposing them allows us to make 
contact with current literature in cosmological perturbation theory^. 

With this decomposition, the Lagrangian can be written as a sum of two uncoupled pieces 
for the spin-0 and spin-1 components. Defining Si = diV, it is 

62C = + £(1) (16) 

where 



£(0) = 2 



(3,S' - mS'd'S^ - P2{d,Sy - mS'd.S' (17) 



is the spin-0 field Lagrangian, and 



£(1) = 2 



piN^ - pidiN^d'Nj (18) 



is the spin-1 field Lagrangian. Note that we have eliminated the cross terms between S** and 
A^*, as well as the (3^ term for the spin-1 component via integration by parts. 



^ We often use the superscripts and to denote spin-0, spin-1 and spin-2 fields and quantities 

associated with them respectively. We have not seen spin-2 fields yet, but we will when we consider 
cosmological perturbations and gravity waves. In cosmological perturbation theory, these are often called 
"scalar- type" , "vector-type" and "tensor-type" perturbations |2l| . To avoid phraseological disasters like 
"scalar perturbations of the vector" , we do not use this terminology though we will helpfully remind the 
reader once in a while. 
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We now reap the first benefits of our diligence by computing the equations of motion. 
Varying the action with respect to S"^ and A^* we get 

- ^^±^^d'S. = (19) 
Pi 

Ni - d'^Ni = (20) 

which are simply wave equations with plane wave solutions 

5'(k) oc e-*^^'"^*+*'^-* (21) 
N\k) oc e-^^*+^''-^. (22) 



Note that the the spin-0 component has a propagation speed of 

(23) 



^(0)2 _ f^l + f^2+ 



Pi 

while the spin-1 component propagates at the speed of light. If we insist that the spin-0 
component does not propagate superluminally, then we obtain the following constraint: 

^1±A±A < 1. (24) 

Since Lorentz invariance is being violated, it is not surprising that the massless spin-0 
field Si may not propagate at the speed of light. Here we have imposed the condition that its 
propagation speed to be less than unity. Since we are already violating Lorentz symmetry, 
one can argue that it could propagate faster than the speed of light. However if we allow for 
superluminal propagation of the fields, in a curved spacetime, the theory a priori does not 
exclude the possibility of the vector fields tilting in such a way that they can fit together to 
form a closed loop, leading to the unphysical flowing of energy around closed timelike curves. 



IV. QUANTUM PERTURBATIONS IN FLAT SPACE 

A cursory look at the Lagrangian Eqn. (fT^ tells us that if (3i < 0, the sign of the kinetic 
energy is wrong which potentially will lead to quantum pathologies. In this section we will 
quantize the fields and show that this is indeed the case, i.e. /3i > if we are to have a 
positive definite Hamiltonian for both spin-0 and spin-1 components. 



8 



The bonus from this exercise is that we obtain a consistent way of choosing initial condi- 
tions for the vector field, a feature we will need when we put this vector field into an inflating 
universe in Section ()VI|) . 



A. Spin-0 

We begin by quantizing the spin-0 component. From the spin-0 Lagrangian Eqn. ()17|) . 

= 2 [(3iS^ - Pid^S^d'Sj - /32{diS'Y - padiS^djS'] (25) 

the canonical conjugate momentum is 

.f'.^ = 4ftS,. (26) 

Promoting the variables into operators Si Si, irf'^ Tff'\ we expand them as a sum of 
ladder operators and their associated mode functions 

^'-Wfnl + "*sSJ'°'i^™.) (27) 

4/3i / ("C^ri-i". + atjS:(»)F:'°>.) . (28) 



-(0) 



The ^g'^Vs are eigenmodes of the Laplace-Beltrami operator d"^ such that d'^Y^\ = —k^Y^\] 
a complete description of these eigenmodes, which are simply representations of SO (3), can 
be found in Appendix 

Following the spin-statistics theorem for integer spin fields, we impose the following equal- 
time commutation relations on the operators: 

[5,(x),5^(x')] = (29) 
[7rf)(x),FW(x')] = (30) 
[^,(x),7r^(°)(x')] = t6{5l-5l')5i. (31) 

Plugging the expansion Eqn. (|27p into the equation of motion Eqn. (|19|). we obtain the 
solution 

4") = Afe-^''' + A^^^e^-^^'^ (32) 
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where k = |k|. If we choose the vacuum state |0) to be the one annihilated by the lowering 
operator ctj^lO) = 0, then we call the term the positive frequency mode and set A_ = 0. 
In order to normalize the ladder operator commutation relations to unity, we choose the 
normalization 



A 



(0) 



1 



Using Eqn. (j27p and Eqn. (j28p . we find that the ladder operators obey the commutation 
relations 

[ag,atg,] = sgn(A)5(k-k') (34) 

[^k'^k'] = (35) 
[atg, atg,] = (36) 

confirming our previous suspicion that (3i > for the field to have the right spin-statistics 
relation. In other words, we need (3i > for spin-0 states to have positive norm (i.e. non- 
ghostlike). Since 4°)^ = (/^^ + p2 + p3)/Pi > for the field to be non-tachyonic and stable, 
we also find that Pi + P2 + Ps > 0. 

We can check that our choice of sign for j3i is the right one by computing the Hamiltonian 



H 



1 



which clearly shows that /?i > for it to be positive definite. 



sgn(/3i)- I dfkicf^k) a>k + « -k«-k + <^(0) 



(37) 
(38) 



B. Spin-1 



The quantization of the spin-1 component proceeds analogously to the spin-0 component; 
there are just more indices to keep track of. As usual, we begin with the Lagrangian Eqn. 

m 



£(1) = 2 

which has the conjugate momentum 



ON' 



(39) 



(40) 
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Promoting the variables into operators Ni Ni, vrf^'' — > vrf^'*, we expand them as a sum of 



ladder operators and their associated mode functions 

^ / d'k (aj;vl"'y<"'. + atsiV;<«)F^'"'.) (41) 
J, / A- (ajiVr"yf + atciV,-l")r:'">,) , (42) 



(27r)3/2 



where are the analogous Laplace-Beltrami eigenmodes for the spin-1 field. The ladder 

operators here are not the same operators from the spin-0 case; they act on a different 
Hilbert space. We have dropped the superscripts on the ladder operators for the sake of 
simplicty. The right hand side is summed over both ±1 modes. As spin-1 fields, we impose 
the equal time commutation relations 

[7V,(x),iV^(x')] = (43) 
[7rf)(x),7r^-(i)(x')] = (44) 
[iV,(x),7r^'(i)(x')] =25(x-x')5^' (45) 

on the operators. 

Inserting the expansion Eqn. 1)4111 into the equation of motion Eqn. ^H^ . and choosing 
the vacuum such that ag|0) = 0, we obtain the solution for the mode function iV^^^'' 

iVf ^) = -^e'^'\ (46) 



noting that the propagation speed c^^-*^ = 1 for the spin-1 field. Here we have dropped the 
negative frequency term, and again normalizing it such that the ladder operator commutation 
relations are normalized to unity, taking into account the fact that we sum over both ±1 
modes in (gT)) and (H^ . 

Using Eqn. (jUJ and Eqn. (jUl), we find that the ladder operators obey 

[ag,atg,] = sgn(/5i)5(k-k') (47) 
[ag, a^,] = (48) 
[atg, atg,] = (49) 
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so that, as in the spin-0 case, /?! > for the spin-1 field to have the right spin-statistics 
relation. 

Finally, we can compute the Hamiltonian 



H 



(fx 



sgn(/5i 



£(1) 



d^kk 



aW + a"^ ca r + 5(0) 



(50) 
(51) 



to show that (3i > for it to be positive definite. 

Thus for the perturbations of the LV vector field to have a consistent, non-tachyonic 
quantum field theory, the values of /3j are constrained; namely /3i > 0, /5i + + /^s > and 
{f3i + + /^s)//?! < 1- There is an additional constraint on /?j that we will find when we 
consider the propagation of gravity waves. We will summarize all the constraints on f3i in 
Section (jVTD|l below. 



V. LV VECTOR FIELD PERTURBATIONS IN AN FRW UNIVERSE 

In this section, we compute the equations of motion for the evolution of the vector field 
in a Friedmann- Robert son- Walker universe, adding in the effect of gravity, which we have 
neglected in the last two sections. The reader may want to skip this section, and refer back 
to it for the equations. 

Since the vector field itself does not contribute any stress-energy tensor without the pres- 



ence of other matter fields in the background 



3, 



it is inevitable that we will have to include 



the perturbation variables of these matter fields. 

The exception here is that of an FRW universe with the vector field and a cosmological 
constant. A cosmological constant contributes energy density to the stress-energy tensor in 
the background which the vector field will track, but it has no perturbations of its own, 
allowing us to solve the vector-gravity theory analytically in Section (jVIj) . 
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A. Background Results and Perturbation Variables 

In linear perturbation theory, the Einstein equation is spht into background and perturbed 
components, both which are solved separately : 

G^y = SuGSfjiv , 5G^^ = SttGJT^^. (52) 

For simplicity, we work with a spatially flat universe; its linearly perturbed FRW metric 
in conformal time is 

ds"^ = a{riY [-(1 + 2<l>)rfr/2 + Bidridx' + (7^^- + 2^7^^- + 2HTij)dx'dx^] (53) 

where $, Bi and the traceless Htij are the four metric perturbation variables. The 
conformal time 7] runs from —00 to 0. 7^,- is the fiat spatial metric in 3 dimensions. 

We have shown in a previous paper jl^ that the background energy density and pressure of 
the vector field precisely tracks those of other matter fields present in a Friedman-Robertson- 
Walker universe, the net effect being that the gravitational constant G^^ is rescaled as first 



noted by 



22 



Gc = (54) 

l + 87rGa ^ ' 



where we have defined 

a = {f3^ + 3(32 + f33W- (55) 
In the same paper, we have also shown that in the Newtonian limit the vector field similarly 
rescales the Newton's constant by 

where 

6 = -f3im\ (57) 

Note that since our experimental tests of the Newton's constant is done in the Newtonian 
limit. Gat is the value we actually measure. As we will see in Section (|VrD|l below, the 
constraints on the theory parameters /3j mean that 

a > (58) 
6 < 0. (59) 
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In other words, the cosmological Newton's constant Gc is always smaller than the measured 
Newton's constant Gn, except in the case where m = 0. 
The background Friedmann equations are 

n = Pm (60) 

(-) = ^a\pra+P^) (61) 

where Ti = a' /a = aH while primes denote derivatives with respect to conformal time t]. 
and Pra are the density and pressure respectively for all other matter fields present in the 
universe which we model perfect fiuid 

T^v = (Pm + Pm)n^n^ + p^gi^u- (62) 

is a unit timelike vector field representing the fiuid four-velocity. These matter fields can 
be dark matter, radiation, or some other exotic isotropic source. 

Using the background result = [m, 0, 0, 0) and the fact that the fixed norm Eqn. (jS)) 
constrains the total length of the vector {u + v)"^ = —m?, we find that the timelike component 
of the vector perturbation is proportional to $: 

u° = $; 63 

a 

i.e. as in the fiat space no-gravity case in the earlier two sections, we only have three spatial 
dynamical degrees of freedom. For computational convenience, we write them as co-moving 
dimensionless variables 

= -V\ (64) 

a 

We then insist that the spatial indices on V^, like those on the shift perturbation i?*, are to 
be raised and lowered using the spatial metric 

V, = -f,jV^ = V\ B,= ^,,B^ = B\ (65) 

Note that this definition means that the original vector perturbation obeys Vi = ma{Vi — Bi). 
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B. The Bardeen Decomposition 

In an FRW background which is spatially homogenous and isotropic, the perturbation 
variables can be decomposed into spin-0, spin-1 and spin-2 components [2^. If further, 
the background matter fields themselves do not break spatial SO (3) symmetry, the great 
advantage of this decomposition is that each of the components are decoupled from one 
another, thus allowing us to solve for them separately. Fortunately, this is indeed the case 
in our theory, since the background vector field only has timelike components (unlike the 



case m 



24| . for example). 

In this section, we will write down the equations for each component. The equations here 
may seem daunting, but they are straightforward to derive: just insert the perturbed metric 
(I53|) and vector field u'^ into equations (@)) and (jHI) from Section (jlTl) and then keep all terms 
;o linear order. The mode equations for the Einstein tensor can be found in the literature 
21| so we will not list them here. We will put these equations to work only in the following 
sections, so one can skip this section and just refer back to it for the equations. 

Using the eigenmodes described in Appendix (^, we can expand the variables as follows 

$ = ^$r(o) (66) 

k 

^ = (67) 

k 

V' = ^ V^^"''>Y'^^"''^ (68) 

k m=0,l 

= ^(^"V^^^'") (69) 

k m=0,l 

= E E H^^"'^Y'^^^'^^ (70) 

k m=0,l,2 

where we have dropped the implicit subscript k on the right hand side. Plugging these back 
into the equation of motion Eqn. and the stress-energy tensor Eqn. (jHl), we obtain their 
corresponding spin-0, spin-1 and spin-2 components. The equations in real space are listed 
in Appendix 



1. Spin-0 "Scalar-type" perturbations 
The spin-0 component of the equation of motion Q is 



a 



{ - Q^n'^ + 6(32 





$ + a^6X + n ' 







for = 0, and 



a 



a a 



\ a 



(5(0)_|/(0))_^ 



(3i + 



a 



for V = i. 

The components of the stress-energy tensor for the spin-0 fields are 

2 



6T 



.0(0) 



5T 



'O(O) 



\ a 



-Pi 



(^(0)_5(0)) 



a 

2 



5r 



:(o) 



(0) 















7^2-2 













+ a 



<t>-n—<^' 



m"' 3 



-2(-) \iP,+Ps)a- 



y».(o) 
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2. Spin-1 "Vector-type" perturbations 



a 
a 



The spin-1 component of the equation of motion Eqn. is 

(5{±i) _ v(±i)) 

The components of the stress-energy tensor for the spin-1 fields are 



^yO(±l) 

^2.o(±i) 



2(^V 

a ) 



-2^n- + 4 - 



-A 



(±1) 



(^(±1) _ 5(±l)^ 



3 



(76) 



(77) 



(78) 
(79) 



3. Spin-2 "Tensor-type" perturbations 



There is no spin-2 equation of motion for the vector field, since it has at most spin-1 
components. However, its stress-energy is non-zero since it contains curvature terms : 



5T 
6T 



'0(±2) 



(±2) 



6T'l^'^ = 



2(^V 
a J 



(80) 
(81) 



VI. LV VECTOR FIELD PERTURBATIONS IN DE-SITTER SPACE 



The paradigmatic theory of early universe is currently the theory of inflation [25[, which 
postulates that the early universe undergoes an accelerated phase of expansion. One key 
to its success is that it provides a mechanism for consistently choosing initial conditions to 
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produce a scale-invariant spectrum of cosmological perturbations that eventually will grow 
into galaxies and things like human beings 3|. 

This mechanism is based on the idea that quantum fluctuations in a near de-Sitter back- 
ground can be "blown up" in their amplitude such that they become classical perturbations 
after they cross the Hubble radius. Thus, the initial conditions for the perturbations are 
quantum mechanical in nature; one chooses a suitable vacuum for the initially quantum field 
when the physical length of the mode is very small, and then follows the evolution of this 
mode until it freezes out. The classical amplitude of this mode is then given to us by the 
2-pt correlation function of the quantum field 2o| . 

In general, inflation is driven by some dynamical field called the infiaton; in the usual case 
it is a scalar field. The infiaton itself has perturbations, so a coupled vector-infiaton-gravity 
theory is not easily solved. Before we do that in the next section, as a warm up we tackle the 
simpler case of inflation driven by a cosmological constant, which has no perturbations. In 
this case, the vector-gravity perturbation equations are analytically solvable in closed form. 

The solutions to the various equations are straightforward, though it takes some algebraic 
manipulation to write them in a useful form; here we sketch the derivation. One can of course 
skip them and just look right ahead at the solution in Eqns. (jH2I), and (jHE))- We will 
turn the order around this time and compute the solutions for spin-2 perturbations first, 
before finishing with the spin-1 and spin-0 perturbation solutions. 



A. Evolution of spin-2 perturbations in de-Sitter space 

Although the LV vector field has no spin-2 component, as we have shown in Section (j^ 
it has a non-zero stress-energy tensor; its very presence will modify the evolution of gravity 
waves. 

Spin-2 equations are gauge invariant'^. By using the i — j component of the Einstein 

^ Since gauge transformations are generated by difFeomorphisms which in turn are generated by a vector 
field, this means that there are no spin-2 gauge transformations. 
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equation, we obtain the equation of motion for the spin-2 gravity wave modes 

1-7 

= lQnG,m\(3i + (33) (83) 

which is the same equation for propagation of gravity waves in de-Sitter space without the 
vector field except that the propagation speed has been rescaled 

cf)^ = j^. (84) 
1-7 

This gives us an extra constraint that /5i + /?3 < if they are not to propagate faster than 
the speed of hght. We can understand this by reahzing the the vector field has changed the 
permeability of "empty space" . This was first noted by the authors of [l3| in a Minkowski 
space setting, whose results for the propagation speeds we also recover for spin-0 and spin-1 
fields in the following sections (by setting a = constant). Note that their "small Cj" limit 
corresponds to IQ-KG^^jSim? <^ 1 here. 

Since gravity wave perturbations freeze out at the sound horizon c^^^ki] = 1, they will 
freeze out earlier compared to the usual case of ki] = 1, leading to a modification of the 
amplitude of the spectrum. We leave the computation of this spectrum to Section (jVIIj) . 



B. Evolution of spin-1 perturbations in de-Sitter space 

Unlike the spin-2 equations, the spin-1 Einstein equation has two gauge degrees of free- 

I I 

dom, and one can use this freedom to simplify the equations somewhat. Following [16], we 
choose 

H^^^^ = (85) 

i.e. we consider the perturbations on hypersurfaces that are shear-free. 

The spin-1 i — j component of the Einstein equation gives us the constraint between the 
shift B^^^ and the spin-1 perturbation V^^\ 

5(±i)=^y(±i). (86) 
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Inserting this into — i component of the Einstein equation results in 



2—W 



a I a 

~~2 



1 

+ 2 



(/3l-/33) + 



a J 
1-7 



a 
a 



y(±i) + 



which after defining 



becomes 



K±i)// 



+ 2 



1 + /53M 



7 



0. 



Here we have defined the useful background quantity 



A 



2W 



a 
a 



2Pm+Pn 



(87) 



(88) 



(89) 



(90) 



using the Friedman equations (jUUI) and (jHH). 

In de-Sitter space the pressure of the cosmological constant is equal to minus the energy 
density, = — Pm, so that A = 0; thus the last term in Eqn. (jHIIjl vanishes leaving us with 
an oscillatory solution for This means that 

aiVo) 1 



kr] 



(91) 



with the propagation speed 



.{±1)2 



(1 - + 



1-7 



(92) 



The normalization is set by choosing an appropriate early time rjo for the mode's "birth", 
and then matching V^'^^'' to the quantum vacuum state nI"^^^ we previously found in Eqn. 
fl4fi|l of Section pVj) . This matching is reasonable because at very short wavelengths where 
presumably when these modes are "created", they do not see the large-scale curvature of 
space time. In this case, we can make this matching when {a{riQ) / k) / ^ 1 where H^^ 
is the de-Sitter radius (or the scale of inflation in a near de-Sitter cosmology). The vector 
norm m is there because we have factored it out in our change of variables in Eqn. fl64jl . 
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Notice that the propagation speed here is not unity as in the flat space case: the presence 
of gravity has changed the value somewhat (although the condition f3i + < implied 
that c^^^-*^ < 1). This is due to the fact that the vector field is no longer propagating 
through "vacuum" , but instead has to traverse through a "medium" of backreacting metric 
perturbations; one can check this by turning off gravity (setting G^, = 0) and recovering 
c2(±i) = 1. (See also Q.) 

The spin-1 component of the shift perturbation sources both E-type and B-type 

)olarization anisotropies of the CMB, the latter dominating the former by a factor of 6 
3|. Thus the quantity of interest that we want to compute is the spectrum of spin-1 
perturbations, given by the two point correlation function 

(0|i?(±i)(x)5(±i)(f + f)|0) =7' / ^^^^^J-e\V^^')\^ (93) 

J k kr 271^ 

But since oc a and a grows exponentially in de-Sitter space, can only decay. 

Unlike the infiaton, the vector field does not freeze out, leaving us with nothing but a bunch 
of complicated equations. In addition, if the infiaton field is a scalar field which contributes 
no spin-1 perturbations, this result is general: scalar field driven inflation with a fixed-norm 
LV vector field does not generate a spin-1 spectrum of perturbations. 

In fact after the end of inflation, the first term in Eqn. (jHUj) will decay as and in 
matter and radiation dominated eras respectively, leaving our oscillatory solution for ^^"^^^ 
unchanged thus even if these modes are created during the CMB era they will decay away ^. 

Thus we have shown that the vector field cannot source spin-1 type perturbations, and 
therefore will not generate B-type polarization modes in the CMB. 

C. Evolution of spin-0 perturbations in de-Sitter space 

Analogous to the spin-1 equations, the spin-0 equations have two gauge degrees of freedom, 
which we fix by choosing 

= , = (94) 

This is technically only true for the usual CMB fluids such as photons, baryons and dark matter. If we 
have other vector sources, the evolution of V^^^^ will be non-trivial. 
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which is sometimes called the newtonian gauge or longitudinal gauge. 

Using the off-diagonal term of the i — j Einstein equation, we obtain the constraint 

which we can then plug into the — i component to obtain 

a-\ak<!>y = - — — — — —\{-l6nG,m^pi + 2-fy 



-(27 + IQnG.a 
where we have similarly defined 



2 - 167rG,m2/3i P * a 



2H^ - 



a 



(96) 



^(^)=aV^^K (97) 

Inserting this result into the space component of the equation of motion Eqn. (ff^ . together 
with the constraint Eqn. ()95|) . we obtain the following wave equation 



^(0)„ ^ (A +/^2 + /?3) ^^2^(o) ^ ^^^(0) _ 

where 



■^2 -^2 



Z, ^ - (99) 

^ /^i(l-7) . s 

' GJG, - SnGcm^pi ' ^ ' 

Comparing ()98|) to the spin-1 solution Eqn. ()89|1 . we see that except for a different propaga- 
tion speed 

c<»>^ , (A±A±A) ,101) 

and the coefficient in front of the last term on the left hand side, they are identical. One 
can check that the constraints on (3i (to be listed below) means that the propagation speed 
c^g'^ < 1 and that it reduces to the flat space case (/?i + /32 + h)/ I3i when gravity is turned 
off by setting G^ = 0. 

In de-Sitter space Pm = —Pm so that A = 0, the last term in Eqn. (jUH)) vanishes and we 
have the solution 

_ a(nn) 1 
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where again the normahzation is chosen to match the quantum vacuum state Sk (Eqn. (|33|)) 
at some very early "birth" time tjq. 

We have shown that the spin-0 perturbation will decay exponentially in an inflationary 
phase and contribute no spectrum of density perturbations on the sky. However, in the 
radiation and matter dominated eras, the spin-0 equation of motion no longer behaves like 
Eqn. (j98|) since it will couple via gravity to spin-0 perturbations of the other fluids, leading 
to complicated coupled equations which cannot be solved analytically in general. 

D. What are the allowed values for /3j? 

Before we proceed to the next Section, let us summarize our flndings on the constraints 
on the theory parameters Pf. 

1. Subluminal propagation of spin-0 fleld : {Pi + P2 + Ps)/^! < 1- 

2. Positivity of Hamiltonian : Pi > 0. 

3. Non-tachyonic propagation of spin-0 fleld : {Pi + P2 + Pz)/ Pi > 0. 

4. Subluminal propagation of spin-2 fleld : A + /?3 < 0. 

Conditions El and El imply that Pi+ P2 + Pz > 0, conditions [Hand El imply that P2 + PZ < 0, 
while conditions El El and El imply that P2 > 0. These constraints form an inflnite length 
wedge in the parameter space spanned by {Pi,P2,P3) as shown in Fig. ((T)). 

The above conditions imply that the background energy density of the vector fleld in an 
FRW background 

Pu = -3aH^ = -3{Pi + 3P2 + PsWH^ < (103) 

is non-positive since a > 0. As discussed in the a companion paper [3], this is not a 
disaster since the background vector has no dynamical degrees of freedom. This is akin to 
the negative energy density of a negative cosmological constant. As shown in the previous 
section, for suitable choices of the parameters Pi, the quantum fleld theory of the propagating 
perturbative degrees of freedom is perfectly consistent. 
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A 




FIG. 1: Figure showing the allowed values for The arrows indicate extension to infinity. The 
slope of the wedge makes a 45 degrees angle with the /?2 — plane. An example of a valid choice 
of parameters is /3i = 1, /?2 = 2.5, /?3 = —3. 

We note here that the conditions above are valid in the regime where 167rG*/3jm^ <^ 1; we 
can understand this by noting that in this regime the energy density of vector field becomes 
small and thus the geometry of spacetime approaches that of Minkowski where we have 
derived our constraints. While we have shown above that the various propagation speeds 
are still subluminal in a de-Sitter background, there is no guarantee that this would be true 
in general. Indeed outside this regime, the modes may become tachyonic^ 

VII. LV VECTOR FIELDS IN INFLATION 

One can try to understand the previous section by analogy to the exact tracker solution of 
the background vector field: the perturbations need to track something. If the background 
has no matter fields and thus is Minkowski, the stress-energy of the vector field vanishes 
identically. Similarly, in the absence of other matter perturbations; the vector field pertur- 
bations follow the behaviour of metric perturbations, which in de-Sitter spacetime have been 



^ We thank Ted Jacobson for pointing this out to us. 
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shown to exponentially decay 

It then follows that if we drive inflation with an infiaton instead of just a plain cosmological 
constant, the vector field might track the evolution of the infiaton perturbations, preventing 
it from decaying away. In this section, we investigate this possibility. 

The simplest models of inflation are those driven by a slowly-rolling scalar field (see for 
example |28|). Unfortunately, unlike the previous section, the equations become intractable 
to solve in exact closed form. Nevertheless, in Section (jVII A|) . we derive approximate so- 
lutions for the evolution of the metric perturbation $. We find that, like its background 
counterpart, the infiaton perturbation acts as a source for the vector field, preventing it 
from decaying away. 

We end our investigation into observability of LV vector fields by computing both pri- 
mordial spectra for the density perturbation and the gravitational wave spectrum in Section 
()VIIB|) . showing that both spectra are modified from the usual values in potentially de- 
tectable ways. 



A. Approximate solution to the Metric Perturbation $ 

The observed spectrum of temperature anisotropics in the CMB [2^ is given by the two 
point correlation function of the density perturbation 5p/p at the time of recombination. 
The density perturbation itself is sourced by a spectrum of metric perturbation $ modes, 
which remain constant until they enter the Hubble radius, at which point they oscillate to 
generate the observed features. 

Thus it is of interest to us to investigate the evolution and the eventual large scale solution 
to $. In this section, we derive the approximate solutions to $, and then in the next section 
we use these solutions to calculate the spectrum of perturbations generated by inflation with 
the presence of the vector field. 

We note in passing that since the infiaton field is a scalar, it sources no spin-1 type 
perturbations and thus the spin-1 perturbation evolves exactly as in the previous section. 
Since we are considering only spin-0 quantities here, we drop the superscript. 
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Consider a scalar field theory acting as the infiaton 



S^= j d''x^[--{d<PY-V{<P) 



(104) 



The exact form of V^(0) is not important in what we are considering; suffice to say that it is 
almost flat (see for example 0]). For this action, the scalar field equation of motion is 



o dV 

h" + 2n(t)' + a^— = 0, 



(105) 



where the scalar is assumed to be spatially homogenous as required by the symmetry of the 
background FRW metric. 

The background stress-energy tensor of this scalar field has the following components 

1 



V 



t: 



4>3 



2 



(106) 
(107) 



Thus if the potential V is very flat, then the scalar field is slowly-rolling cfp'' <^ Va?' and thus 
Pcj, ~ —P(t> leading to an accelerating phase of cosmic expansion. 

Perturbing this field (j) — > (f) + 50, we find its linear order spin-0 stress-energy tensor is 



Y 



1 - 



'k5(i)Yi 



ST 



<t>3 



kl2 



+ 



dV 



¥6' 



(108) 
(109) 
(110) 



Combining these with the vector stress-energy tensor components Eqns. ()73|) . ()74|) . (f75j) 
from Section ()V B 1|) and the Einstein tensor, we obtain 



dV 



\ a , 



+MGcl3im^ (kV + knV) - aHkV 



{kn^-k<iJ') = ^(l)'5(l)-aAV + Pim^a'\ak^y - Pm^ — 



a 



(111) 
(112) 



AttGc {(I)'6(I)' - ^a^c 



\ a . 



^''^^'"'(A + /?2 + /?3)A:^($ + vi/) 



7 



(113) 
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which are the — 0, — i and i — i components of the Einstein equation respectively. 

Our goal is to write down the two evolution equations for the coupled $ and ^ = aV 
fields, by eliminating and 6(j). We first replace every instance of \1/ with $ and C, by using 
the constraint equation (jH^). The first equation is then simply the equation of motion (f7^ 



e+ 



a(l -7) 



/3im2 + 70; Pim'^ + 70; 

where now the background quantity A is given by 

\ a 



A 



2W 



Piw? + a 
Piw? + 7a 

2P<I>+P<t>' 



A;(a$)', 



;il4) 



(115) 



Next we obtain the second evolution equation by first adding Eqn. pil|l and Eqn. piHj) 
to eliminate 6(f)'. To write 50 in terms of the $ and ^ variables, we use the equation of 
motion Eqn. ((721), combined with Eqn. ()112|) . and the background equation ()1U5|) . to obtain 
the following relation 



- %TxGc-r-a^6(\) 



A 



[H^ + $') - 



2i! 

k a 



87rG,(A + /32 + /33)m2fc-^A)^ 

k J a 



(116) 



Putting everything together, and eliminating the dV /dcf) term in Eqn. (jll6|) using Eqn. 
()105p . we get the second evolution equation 



, a"\ A' 
-2n^ + 2 — - a—n 



A 



+87iGc 

7 



+ 



-- 0. 



7 



k^ + 8nG,m\f3i + ^2 + /^s)^^- 
a 



A a a a 



A a 

(117) 

To further simplify our analysis below, we set A — > 0, i.e. the background is very near 
to de-Sitter space. This means that a"/a ~ l/rj'^, and a'/a ~ A'/A ~ l/?7, with ~ meaning 
"scales as". 



1. Short wavelength solution 



Consider the evolution of the perturbations at short wavelengths krj ^ 1. Unlike the de- 
Sitter case, V = ^/a no longer decays away since the presence of the inflaton perturbation 
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acts as a source. We can see that by using Eqn. (jll2j) and doing similar computations to 
Section ()VI C|) to obtain 

e , (/3i + ^2 + /33) ,2e Am^ + a 

— + ^ - = -5 571 r47rGcfc0 S(f). 118 

a Z2 a pim^(l — 7) 

The source term on the right hand side is proportional to k, and thus the inflaton per- 
turbation only sources the vector perturbation at small scales. In contrast, the de-Sitter 
background case Eqn. is always sourceless. 

To investigate the evolution of $ however, it is more convenient to use Eqn. (|114j) which 
we can solve (recall that we have set A = 0) 

e = [D^ cosicskr]) + D2 sm{6skr])] + f ^ / k{a<^)dr], (119) 

Pim^ + 07 J 

where 

c:- ^ 'fty'+'^'>'"' (120) 

Pim'^ + 07 

and Di and D2 are integration constants. The bracketed first term on the right hand side 
is just the usual oscillating solution (albeit with a modified propagation speed) we obtained 
in Section (jVI Eqn. ()l()2j) . Since V = ^/a, the oscillating terms decay away like their 
de-Sitter space counterpart. Note that this solution is valid for all wavelengths, ii A = 0. 

Plugging solution (jll9p into Eqn. ()117|) . and dropping all terms that scale smaller than 
k"^, we obtain the short wavelength evolution equation for $: 

$" + [ 87rG',m^(A+/?2 + /33)(l + aMm^) ^ 1 ^ ^ 
[ (1-7) J ' 

whose solution is a plane wave with a propagation speed of 

2 _ 87ram^(A + /32 + /?3)(l + «Mm^) , , 

- (1-7) ^ ^ 

Setting = recovers the usual solution for inflation without the presence of the vector 

field. 

As expected, the metric perturbation has its propagation speed modified, consistent with 
the argument that the ever-present vector field has changed the permeability of spacetime. 
Note that the conditions listed in Section ()VID|) mean that the propagation speed is less 
than unity, so $ propagates subluminally. 
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2. Long wavelength solution 



Now let us consider the evolution of long wavelength kr] <^ 1 perturbations. In this 
regime, Eqn. ()114j) becomes 

r = 0, (123) 



which has the simple solution 



(124) 



where and D4 are the constants of integration for the decaying and the constant terms 
respectively (recall that i] runs from —00 to 0). Plugging this back into Eqn. ()117|) we get 



A'' 




m - A 







-27V + 2 



a 
a 



A' 
~A 



$ = 0, 



(125) 



which has the following solution in an inflationary background 



adt 



(126) 



where the cosmic time J dt = J adrj. At the end of inflation when the universe enters a 
radiation dominated era, a oc t^''^, so (1 — Ha~^ J adt) ^ 2/3, thus the frozen-out spectrum 
of $ is then 



(127) 



The value of depends on the perturbation sources at the same freeze-out time, i.e. at 
the moment the mode leaves the Hubble radius^. We will compute this value in the section 
below when we derive the perturbation spectra. 

One important result is that, although we have an extra field on top of the infiaton, no 
isocurvature perturbations are produced at superhorizon scales, since Eqn. (jl25|) has no 
source terms on the right hand side j^]. This is an interesting departure from usual multi^- 
scalar field models of inflation, where isocurvature perturbations are a generic feature j3lj 
(see also 1321). 



^ Technically speaking the mode freezes out when it leaves its sound horizon, i.e. at Cgkg = 1. 
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B. Perturbation Spectra and the Consistency Relationship 



Observables on the CMB are ultimately sourced by the frozen-out spectra of super-Hubble 
perturbation modes reentering into the Hubble radius after their sojourn outside the Hubble 
radius. In this Section, we compute the perturbation power spectra for both the spin-0 
"scalar" and spin-2 "tensor" fields. 

We first tackle the spin-0 "scalar" perturbation spectrum. From Eqns. ()112j) and ()72j) . we 
replace with Eqn. pi8j) . and then using the non-decaying part of the short wavelength 
solution for Eqn. ()119|) . we get 



(a$)' = AttG 



1 + 7a/ {Piw? 



where 



1-7 

(/?l+/52 + /?3K 



(128) 



(129) 



/5im2(l — 7) 

Eqn. p28|l is the solution for 50; we solve the equation of motion for $, and then insert 
the result into Eqn. ()128j) to find 6(f). In the long wavelength regime fcr] -C which we are 
interested in, the term drops out. 

We want to match the long wavelength solution of Eqn. 
the long wavelength solution of $. To do that, we multiply Eqn 
time derivative of the result to get 



to Eqn. ()126|) . which is 
with a and take the 



a J 



^2 H 



(130) 



where we have used the background relation —T-C'+Ti^ = 47rG'c0'^ and / a'^drj = J adt ^ a/H. 
This means that Eqn. ()12(j|l becomes 



a 



which we can match to Eqn. 



D 



to obtain the value for 

1 + 7a/ {Pim?) a5(j) 

' ~ 1-7 0^' 



(131) 



(132) 



Thus at the end of inflation, using Eqn. ()127j) . the metric perturbation $ has the following 
value: 

'1 + 7a/(/5im2) 



3 6' 



7 



(133) 
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Now using the usual result for the power spectrum of the inflaton 50, Vs,j, = (50^) = H"^ / {2k^) 
at large wavelengths, we obtain the spin-0 "scalar" power spectrum 



SttG, 1 



9 e 

Here the inflaton slow-roll parameter^ 



1 + 7a/(/3im2) 



7 



2 



(134) 



,.GA^ . (135) 



Comparing this to the case without the vector field (see for example |33|), the amplitude of 
the perturbations has been rescaled by a factor of 



7^$ G* 



1 + 7a/ (/5im^ 



(136) 
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where P$ is the power spectrum without the vector field. One can check that by setting 
= 0, we recover the usual relation for the power spectrum p34j) where P^/P^ = 1. Note 
that the index of the power spectrum does not change since vector field background precisely 
tracks that of the inflaton. 

As an aside, we comment on the use of the relation 6(f) ^ H above. One can argue 
that there is no reason to believe that the inflaton perturbation at small scales will evolve 
similarly to the usual no- vector-field case. However, since there is no direct coupling between 
the inflaton and the vector field, any effects that the inflaton perturbation will experience 
due to the vector field has to be through the metric perturbation which is at most a second 
order effect. 

This rescaling of the spin-0 amplitude is degenerate with the slow-roll parameter e, thus 
from the standpoint of the spin-0 perturbation alone we cannot disentangle them. Never- 
theless, there is hope for us to distinguish between them if we can observe the gravitational 
wave spectrum, as we will see later. 



From Eqn. (|82p and following 
'tensor" perturbations to be 



we can then compute the power spectrum for spin-2 



A;3 



^ We have used the potential- independent definition of the slow-roll parameter. See |30| for a discussion on 
the difiference between this parameter and the V dependent definition e ~ (WttG^,)'^^ {{dV /d(f>)/V)^ . 
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Since for non-super luminal propagation of the gravity waves we need 7 < 0, the spin-2 
"tensor" power spectrum is boosted, making the tensor spectrum possibly more accessible 
to observations (currently we have no detection of a spin-2 spectrum js^). Of course, this 
value is degenerate with the Hubble parameter H so measurement of the spin-2 spectrum 
does not allow us to put a constraint on the LV vector field parameters. 

If we assume natural C(— 1) values for /3i-|-/?3, then 7 ^ {m/mpY where = l/(87rG'jv). 
Thus if m is of order rrtp which is to say that Lorentz- violating effects are a symptom of new 
Planck scale physics, then there will be a boost in our chances at observability of gravity 
waves in the CMB. However, as shown in primordial nucleosynthesis constraints m to 
be an order of magnitude smaller than rrtp. 

Before we close our investigation of LV vector fields in an inflationary universe, we note 
that the presence of this vector field has violated the so-called inflationary consistency re- 
lation. One can show that any single slow-roll scalar field model of inflation with slow-roll 
parameter e, the spectral index for the spin-2 spectrum ut = —It. But the ratio of the 
amplitudes for the spin-2 and spin-0 spectra is proportional to e, so this relation states that 

^^-(2) 9 

nr. (138) 



P$ 2 

But here, with the vector field permeating the universe, from Eqns. fll34|) and p37j) we 
obtain the modified relation 

!Vi ^ _£ (l + 8xG.)(l-7P/' 
■P* 2 |l+7a/(Am2)12 

In other words, if we believe in the single slow-roll scalar field model of inflation, then 
measuring ut will allow us to put a constraint on the parameters of the theory. 



VIII. CONCLUSIONS 



We have shown that perturbations of a fixed-norm, timelike vector field have a consistent 
quantum theory. If we insist on positivity of the Hamiltonian and non-tachyonic propagation 
of the quantized field, then we find that /5i > for the former and /?i + /52 + /^a > for 
the latter. In addition, non-superluminal propagation of the spin-0 field and gravity waves 
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requires that {Pi + (32 + Pz)/ Pi < 1 and /3i + /^a < respectively. These facts imply that the 
sign of the background energy density of the vector field is non-positive. 

We then considered the evolution of the vector perturbations in a FRW universe. Spe- 
cializing to a de- Sitter spacetime, we find that the perturbations decay away exponentially. 
We argued that this is because there is no other matter perturbations for the vector field to 
track. 

Following that, wc investigated its evolution in the presence of an inflaton field. As 
expected, the vector perturbation is now sourced by the inflaton perturbation, making it 
non-decaying and thus imparting non-trivial effects on the metric perturbations. We derived 
the generated primordial spectrum for both spin-0 (approximately) and spin-2 (exactly) 
perturbations, and showed that their amplitude is rescaled. In particular, we noted that 
if the norm m of the vector field is large enough, it might lead to a boosted spectrum of 
B-mode polarization in CMB. We also showed that the presence of the vector perturbations 
also violates the infiationary consistency relation. 

Finally, the ultimate test of the observabihty of LV vector field induced modifications in 
the CMB lies in our ability to compute the evolution of perturbations from Hubble reentry 
through to the surface of last scattering. This entails solving the multicomponent Boltzmann 
evolution equations, adding the vector perturbation on top of the usual photon-baryon-dark 
matter fluid, an endeavour that can only be done numerically. 
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APPENDIX A: EIGENMODE DECOMPOSITION 

In this paper, we have chosen to work in eigenmodes of the irreducible representations of 
SO (3). This formahsm was first apphed to cosmo logical perturbations by Bardeen [2^. We 
present a short review in this appendix for completeness and also to pin down conventions, 



2l|. 



specializing to a spatially flat universe. For a good comprehensive review, see 

In a background FRW spacetime that is foliated into spatial hypersurfaces labeled by 
cosmic time, the linear perturbations are just functions living in a spatially S0(3) and trans- 
lationally invariant background. If further, there are no background quantities that violate 
this symmetry (e.g. a background vector field with spatial components or an inhomogenous 
scalar field), the perturbation variables decomposed using these eigenmodes are decoupled 
from each other. 

The spin-0, spin-1 and spin-2 orthonormal modes are eigenmodes of the Laplace-Beltrami 
operator 

^2^(0) ^ ^2y(0) ^ Q (Al) 

d^Y}^'^ + erl^'^ = (A2) 
d'Y^'^ + kXf'^ = 0. (A3) 

These modes obey the following conditions 

5.yi(±i) = , diV'^^^^^ , _ (A4) 

i.e. spin-1 modes are divergenceless, and the spin-2 modes are traceless and transverse. With 
these conditions, we can derive their explicit form 

r.(°)(x) = e**^-^ (A5) 

Jjp 

ri(±i)(x) = -^{kinp±ikpni)e'''-'' (A6) 

Kj^(±2)(x) = -Jh-'[e''P{kinp±zkpni)][e^'~'{krn,±iksnr)]e'^^^ (AT) 



where n is a unit vector such that k ■ n = 0, e^^'^ is the Levi-Civita symbol and /c = |k|. Here 
we have restored the k subscript and the function argument (x) to the mode function. 
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Using these basis modes, we construct other modes. Following the convention of j21| . 

F.W = -A;-i9,F(°) (A8) 
= (k-'d.d, + i7.) Y^'^ (A9) 
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In particular the following identities are especially useful when deriving the commutation 
relations and the Hamiltonian in Section (IIVII 



F.f = -Y^% (All) 
y.f ^ = -Y^^^- (A12) 

We end this appendix with a warning that different authors often use different conventions 
for the explicit mode functions. In this work, the explicit mode functions never appear. 



APPENDIX B: EQUATIONS OF MOTION IN REAL SPACE 



In this appendix, we collect the real space equations for reference. 
The perturbed component of the equation of motion Eqn. Q is 



-6(/5i + 2/52 + /33)H'$ + 6/^2 



+H 



(2A + P2)^^V' + + d,B') + 3/^2$' + 3(2/3i + /?2 + 2/^3)^' 



(Bl) 



for the u = component and 



a a 
m? \ a 



a 
a 



a 



(/?i + ^)9.$ + 2/?i(V/-50 



+ 2(-/53 + f3i)d[jd^B,] - (3,d,dW, - {(32 + (^3)did^V^ 



a 



-(A + (33)^,h'^^ + (3id.<l>' - —d,^' - (3i{B'/ - vn = 



(B2) 



for the u = i component. 
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The stress-energy tensor has the following components: 



ST', 



-3(/3i + 3/32 + /33)^'$ + Pia~^ \adi{B' - V') 



+ (/?! + 3/^2 + P3)n{diV' + 3*0 + (3,d^^} 



m 



-2(/3i + 3/32 + /33)^' + (3/32 + Pz) 



-/3ia-2 [a^Vl - B[)\ - (5^a-\adi^)' 
+^(-/3i + (^z) [d\B, - V,) - d,&{B^ - V^] 



m 



sr. = 2— (/3i + 3/32 + /33) 



\ a , 



m - iPi + 3/32 + P3)ws: 



j 



+a 



(B3) 



(B4) 



(B5) 
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